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OPTIMAL  DESIGN  OF  PLL  WITH  TWO  SEPARATE  PHASE  DETECTORS 


A.  Heiman  and  Y.  Bar-Ness 


ABSTRACT : In  order  to  improve  the  operation  of  phase  locked 

loops  in  both  the  acquisition  and  tracking  modes,  loops  with  two 
phase  detectors  are  considered. 

The  optimization  problem  of  choosing  the  phase  detectors  in 
order  to  minimize  the  variance  of  loop  phase  error  is  formulated. 
The  optimization  problem  is  solved  under  several  approximations. 
Despite  the  approximation,  the  resulting  nonlinearities  can  be 
implemented  and  results  in  a PLL  with  improved  performance  in  both 
acquisition  and  tracking. 
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1.  Introduction. 


In  designing  phase  locked  loop  system  (PLL)  one  is  usually 
faced  with  different  requirements  such  as  minimum  phase  error 
variance,  small  acquisition  time,  maximum  locking  range,  etc. 

With  the  usual  way  of  designing  these  systems  one  cannot  adequately 
cope  with  these  different  requirements  and  the  total  performance 
often  is  not  good  enough. 

In  this  paper,  we  will  show  that  by  adding  another  phase 
detector  the  system  can  be  provided  with  enough  additional  freedom 
to  be  able  to  better  deal  with  these  conflicting  requirements  and 
yield  smaller  errors.  In  order  to  find  good  forms  for  the  phase 
detectors  we  will  use  optimization  theory.  The  detector  function 
which  we  obtain  might  not  be  strictly  optimum,  owing  to  several 
approximations  in  the  analysis  which  we  are  forced  to  make.  Yet, 
as  is  clearly  shown  by  our  simulation,  the  resulting  system  is 
superior  to  regular  systems  in  current  use.  This  study  is  believed 
to  be  the  first  which  exploit  optimization  theory  in  this  way. 

Owing  to  the  use  of  additional  phase  detector,  the  problem 
under  consideration  is  one  of  (near)  optimality  designing  a PLL  with 
two  phase  detectors.  The  equivalent  model  of  such  PLL  is  given  in 
Figure  1 and  can  be  considered  as  a regular  PLL  (see  Figure  2) 
with  H(p)g(<j>)  being  replaced  by  the  sum  Hj(p)f($)  + H2  (p)  f j (♦ ) , 
where  Hj(p)  and  H2(p)  are  linear  filters.  Such  a split  of  the 
phase  detector  and  the  loop  filter  is  shown  to  be  useful  in 
designing  the  system  to  perform  well  in  both  its  acquisition  and 
tracking  modes  and  hence  has  a superior  performance. 

Consequently,  the  design  goal  of  the  paper  is  to  find  the 
phase  detector  functions  fj($)  and  f(4>)  that  (approximately) 
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minimizes  the  variance  of  the  phase  tracking  error.  Since  the 

analysis  of  this  problem  with  general  linear  filters  Hj(p)  and 

H.(p),  is  quite  complicated,  only  the  second  order  case  is 

1_H0 

considered  in  this  work:  i.e.,  (p)  = and  I^fp ) = - . 

Clearly  when  fj(4>)  = f ) the  system  reduces  to  a standard 
second  order  PLL.  Hence,  with  two  phase  detectors,  we  obviously 
have  more  freedom  in  the  design.  It  is  shown  that  our  solution 
involves  the  use  of  an  extra  design  parameter  (termed  Tj).  This 
parameter  enables  the  design  of  a PLL  having  minimum  variance  as 
well  as  locking  time  smaller  than  some  prescribed  value. 

In  Section  2,  the  probability  density  function  (p.d.f)  p(4>) 

of  the  phase  tracking  error  is  derived  as  the  solution  of  a 
Fokker  Planck  equation.  Simplified  form  of  this  function  is  obtained 
using  some  approximation  and  assumption.  In  Section  3,  the 
minimization  problem  is  stated  together  with  its  appropriate 
constraints  and  augmented  performance  index.  The  (approximately) 
optimal  phase  detector  functions  which  are  the  solution  of  the 
related  Euler  equation  are  presented  in  Section  4,  together  with 
the  corresponding  optimal  p.d.f.  Section  5 includes  some  remarks 
on  the  optimal  solution;  in  particular,  the  effect  of  the  design 
parameter  T^  on  the  system  locking  time.  In  Section  6,  the  optimal 
phase  detector  function  is  firstly  calculated  and  graphically 
presented  for  Tj  = 0.  Then  for  T^  t 0,  we  compute  and  present 
graphically  the  dependence  of  c$lmin  on  the  input  signal  to  noise 
ratio  (SNR)  and  on  the  design  parameter  T^ . Using  simulation 
results,  the  actual  dependence  of  the  locking  time  on  this 
parameter  is  also  presented.  This  enables  us  to  suggest  a way  of 
designing  an  optimal  system  having  a locking  time  smaller  than 
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certain  prescribed  value. 

Since  the  optimal  solution  depends  on  the  frequency  error 
between  the  incoming  signal  and  the  VCO  output,  it  is  plausible 
to  check  the  sensitivity  of  the  optimal  system  to  a change  in  this 
frequency  error.  This  is  done  by  simulations  from  which  it  is 
concluded  that,  with  the  frequency  error  unknown,  it  is  best  to 
design  for  the  maximum  of  its  expected  value. 

Finally,  in  Section  7 practical  implementation  of  the 
phase  detector  function  is  suggested,  using  the  first  term  of  the 
Fourier  series  expansion  of  the  actual  optimal  functions. 
Simulation  results  with  the  "sub-optimal"  phase  detector  functions 
justify  the  usefulness  of  proposed  approximation. 


2.  System  Equation  and  Tracking  Error  Probability  Density  Function 
Consider  the  problem  of  selecting  two  phase  detector  func- 
tions, fjC^)  and  f($),  such  that  the  variance  of  the  phase 
tracking  error 

al  = f *!(♦), *(♦))<«  CD 


is  minimized  where  p (4> , f j (<t>) , f O) ) is  the  probability  density 
function  (p.d.f)  of  the  error  (modulo  2it)  and  <p  is  the  expected 
value  of  <j> . 

Following  Lindsey  [1]  and  Viterbi  ( 2 ] in  assuming  that  the 
input  filter  bandwidth  is  much  wider  than  the  PLL  bandwidth,  we 
have  that  N^ft)  and  (t)  of  Figure  1 are  zero  mean  white 


Gaussian  noises  of  spectral  densities  (one-sided)  Nn  and  N. 


respectively.  This  assum 


inti  on  enables 

*8  12 


the  usage  of  the 


11  1U) 
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Fokker-Planck  equation  from  whose  solution  the  p.d.f  of  the  error 
p (<}> , f i (<J> ) , f (<t> ) ) , can  be  obtained. 

It  is  easy  to  deduce  from  Figure  1 that  the  system  equation 
is  given  by 

3T  ' fl0  ‘ H0(AKf(*>  * KNf)  - ITT^p  * KNfj)  <2) 

where  we  assumed  that  the  input  signal  includes  phase  modulation 
0 C t ) = Aflt  + 0Q,  Aq  is  constant  and  A is  the  signal's  amplitude. 
Also,  K = KyB  where  B is  the  amplitude  of  the  voltage  controlled 
oscillator  (VCO)  output  signal  and  Ky  is  its  gain  in  radian/ sec/ volt . 
The  loop  filters  are  Hj(p)  = HQ  and  ^(p)  = (1  -Hq)  / (1+tjp) , where 
p = d/dt  is  the  Heaviside  operator. 

To  write  the  system  equation  in  state  space  form,  define 

yo  ‘ yi  • - nTTT  tAICW  * ™gl 

and  hence 


equation  with  respect  to  y^ 
density , 
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yields  equation  (5)  for  the  marginal 


8(r(y0)p(y0)) 


3 2p  Cy0) 


9y 


T 


= 0 


-n  < y0  < 


(5) 


where , 


r(y0)  • -4(A0-H0AKf(y0)  + E(yi  | yQ) ) /H2K2Ng  . (6) 


The  solution  of  (5)  can  be  written  as 


P(y0)  = c0exp  [ ' 


(l  + C, 


z 

r (x)dx] dz) 


(7) 


where  c^  can  be  obtained  from  the  boundary  condition  p(ir)  = 
p(-n)  and  cQ  is  a normalization  factor.  Notice  that  p(yQ) 
depends  on  the  conditional  expectation  ^Cy^ly^)  which  requires 
for  its  calculation  the  knowledge  of  pCyg.y-^).  In  order  to  proceed 
from  this  point  it  is  necessary  to  make  some  approximation  and  we 
will  follow  the  approach  used  by  Lindsey  which  approximated  ECy^yg) 
by  the  right  hand  side  of  equation  (8)  and  write 

AK(1 -Hg) (f x (yQ)  - 71)S{  (0) 

ECyjIyg)  - 2 i AKCI-HqJT.  (8) 

2Tl°f1 

2 

where  7^  and  are  the  mean  and  variance  respectively  of  the 

random  variable  f,(<fr),  and  Sf  (w)  is  the  spectral  density  of  the 

1 *1 

random  process  ( (♦ ) - 7j).  Using  equation  (8)  in  (6)  we  obtain 
after  some  algebraic  manipulations 
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r (x) dx  = - aF(<|>)  + 3F C<1> ) / crf  - 


(9) 


where , 


dv  dFi 

§ - ar  - 

a = 4A/ll0KN,  6 = 2A  (1  -HQ)  Sf  (0)  / T jNqHJk 


and 


81  ' |fl0  • AK(1-H0)TjW)(1  . eH0/«ci-H0)T|i)i . (10) 

In  order  to  facilitate  a solution  we  need  to  make  some 

additional  assumption.  The  restrictions  to  be  imposed  by  these 

assumptions  will  make  physical  sense.  In  fact,  we  will  assume 

r<t> 

that  f.  (♦)  and  f(<f>)  are  such  that  r(x)dx  = 0 for  all 
1 1 -<P 

ip  e [-7i, it].  This  implies  that  c^  = 0 and  p(<J>)  is  symmetric 
and  given  by 

r<P 

p(<l>)  - cQexp  [ - J r (x) dx]  . (11) 

Furthermore,  to  satisfy  the  condition  on  r(x),  we  will 
assume  (see  equation  (9))  that 


«F(*)  - BFj (♦)/<*£  - «F(-*) 


BF1(-#)/«£  . 


C12) 
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One  way  of  satisfying  equation  (12) 
even  functions.  But,  then,  f(9)  = 


is  to  make  F(9)  and 
and  fjC*)  = - 


Fj  (♦) 
dFjf*) 

3? 


are  odd  and  since  p(9)  is  symmetric,  we  have  7^  = 0.  Hence, 

equation  (10)  yields  Bj  = aA^/AKHp  and,  therefore,  with  Bj  = 0 

it  is  necessary  to  have  A^  = 0.  Since  we  do  not  intend  to  restrict 

our  problem  to  AQ  = 0,  we  will  choose  F(9)  and  F ^ (4> ) such 

2 

that  their  weighted  sum  aF(9)  - 6F.  (<(>)/of  is  an  even  function 

1 *1 

and  that  Bj  = 0.  It  is  shown  in  the  sequel  that  this  enables 
us  to  treat  the  case  of  A^  f 0. 

It  is  important  to  notice  that  in  contrast  to  the  standard 
PLL  case  where  the  one  phase  detector  function  is  usually  chosen 
as  an  odd  function,  in  order  to  make  p (<p ) symmetric  [1]  , [3]  , [4]  , 
for  the  two  phase  detector  case  F ' (4> ) and  F ^ (4> ) , are  not  odd. 

In  fact,  when  Afl  f 0,  p(9)  is  not  symmetric  with  respect  to  the 
steady-state  error  9 7*  0 for  the  standard  PLL  even  though  the 
phase  detector  is  an  odd  function.  In  view  of  this,  the  weighted 
sum  of  the  two  detector  functions  is  chosen  to  be  an  odd  function 
so  as  to  ensure  that  p(9)  is  symmetric  even  when  AQ  / 0. 

Finally,  when  equation  (12)  is  satisfied, 

aFOn-FjC-n/Of 

P(<1>)  a jj*  e 1 (13) 


where  D*  is  a normalization  constant.  That  is,  p (4> ) depends  on 
2 

of  , which  itself  depends  on  p(^)  by  its  definition.  Hence,  in 
*1 

order  to  obtain  p(4>)  be  well  defined,  we  must  first  solve  the 

2 

following  simultaneous  equations  for  a, 

*1 


8 


it  aF(«)-BF1(«)/o2 

rf  = I * 7i)2e  Jd<<) 

1 J - 71  1 1 (14) 

TT  «F(4»)-  F^/O* 

^1  = Jpf  J fl  ^)e  ^d$  • 


2 

The  fact  that  p (4> ) depends  on  of  complicates  the 

2 1 

search  of  the  minimum  of  o^  which  will  be  discussed  in  the  next 

section,  furthermore,  8 depends  on  Sf  (0)  and  since  the  latter 

rl 

cannot  be  explicitly  calculated,  a heuristic  approach  (see  Section  5) 

is  used  to  argue  that  Sf  (0)  is  almost  independent  of  f,  (<j>) 

rl  1 

and  hence  that  8 is  approximately  a constant.  This  is  supported 
by  simulation  results. 


3.  The  Minimization  Problem 

The  problem  is  to  find  f(<p)  and  fj(<f>)  so  as  to  have 

minimum  value  for  of 

<P 

oiF(<f>)  - 8F,  (4>)  /of 

7 1 f ? 1 t . 

% = jpr  J ^ e 1d<P  (15) 

dF  dFl 

where,  f(4>)  = * 3^  » f j (♦ ) a - 3^—  . Since  p(4>)  in  equation  (13) 
was  obtained  using  the  conditions  of  equation  (12)  the  optimal 
problem  must  satisfy  the  constraints 


H 3 

[Aq  - AK(1-H0)T.(1  ♦ ^ J- 

0 0 0 1 a (1-Hq) Oj 


■)] 


(16) 
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«F(4>)  - BFj(^)/a^  = F C _ ) 


8F.  -n  < <p  < it 

lj 


(17) 


1 

D* 


7,  «F(4>)-F2(4>)/a2 

e 1 

. -TT 


dy 


1. 


(18) 


2 

The  value  of  of  is  obtained  from  equation  (14). 
rl 

It  is  possible  to  show,  using  equations  (16),  (12)  and  the 
symmetry  of  p(4>)  that 


a-typ;  * h0f’  • - ^ 


(19) 


This  requirement  on  the  phase  detector  functions  will  prove  to  be 
of  great  importance  in  the  system  acquisition  mode.  In  fact, 
t >>  one  may  deduce  from  equation  (3)  that  FJ(4>)  is  stationary 
and  FJ  is  constant  and 


- Xj  = (l-H0)AKF|(<f).  (20) 

Now  averaging  equation  (4a)  and  using  = yQ,  we  get 

3F  * Ao  * HoAKF'(*’  * *i-  12U 

Substituting  equation  (20)  in  equation  (21)  and  using  equation  (19) 
we  end  up  having  ||  -►  0 regardless  of  the  value  of  AQ . This 
means  that  by  satisfying  condition  (19)  the  system  ends  up  locking  for 
any  Ap  and  the  acquisition  range  of  this  general  PLL  is  infinite. 
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This  is  certainly  not  the  case  for  the  standard  second  order 
PLL  where  the  acquisition  range  is  limited  and  depends  in  [1] . 

To  conclude  this  discussion,  we  notice  that  equation  (19) 
can  be  written  as 


0 


jpr  j ^ [HQF  ' (<t> ) + (1  *H0)  (*)J  exp  (aF(*)  - BFj  (< P)/a2f  )d4>  (22) 


which  is  equivalent  to  equation  (16)  . This  together  with  equations 
(17)  and  (18)  form  the  constraints  of  the  optimization  problem. 

Using  Lagrange  multipliers,  the  augmented  performance  index 
associated  with  this  problem  is  given  by 


A ctF(v)-  F 1(^)/a‘ 

14>2  + AjO^F'f*)  + (1-H0)F'(4>)  + j^)Je  1 


J = 


(23) 


-it 


exp  [otF  (<)>)  - 0F.  ) d4> 

1 rl 


• where 


is  the  multiplier  and  the  value  of 


is  given  by 


Fj  ) 2exp  [cx F (4> ) 


]d* 

1 *1 


(24a) 


F{  - - , where  a ■ (1-Hq)  + — . (24b) 

fl 

Although  the  constraint  (17)  has  not  been  explicitly  considered  in 
deriving  equation  (23),  it  can  be  shown  that  it  is  satisfied  in  the 
optimal  solution. 

Because  of  shortage  of  space,  we  will  not  write  the  Euler 


11 


equaton  which  F(<t>)  and  F ^ ) must  satisfy,  (for  full  detail 

see  [3]),  and  we  only  mention  the  boundary  condition  that  these 
functions  must  satisfy: 

Fj  (TT)  - FjC-tt)  £ Tr  F'(TT)  = F'C-ir).  (25)  . 

For  standard  PLL,  = 0,  while  with  the  two  phase  detector  PLL, 
we  can  choose  arbitrarily.  This  adds  a degree  of  freedom  which 

will  prove  to  have  an  important  effect  on  the  locking  time.  Therefore, 
we  can  design  the  system  for  minimum  variance  and  simultaneously 

put  some  requirements  on  the  locking  time. 

4,  The  Optimal  Phase  Detector  Function 

The  system  Euler  equations  are  second  order  differential 
equations.  However,  since  we  are  interested  only  in  F'(<£)  and 
FJ(4>),  we  need  to  solve  first  order  equations  only.  From  this  we 
can  obtain  the  two  optimal  phase  detector  functions 

F1CW  ' Fi  = IT^tT  V -yAKa")  (a"-14T)  1-1  + ^ + ,r2/sinh2xJ  (26) 

and 


F ' (<j> ) 


a-l+H 


ITT  FlC4>)  72'  _2 V ."2 

0 +7r  /sinh  x 


(27) 


which  depend  on  the  system  parameters  as  well  as  on  the  parameters 
a and  x.  These  parameters  are  the  simultaneous  solution  of  the 
two  nonlinear  equations, 
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a V2"  A0 

7 (V2  + A0/AKa^  " Ca-1+H0)  + 


(28) 


-2(1  + 


(V2*  + A0/AKa)^(a-l+H0)«  sinh2  sinh2x 


BH, 


) 


TT 


+1=0.  (29) 


The  phase  detector  functions  are  (approximately)  optimal  in  the 

2 j i 

sense  of  minimizing  and  ^ ■+  0 for  every  A^. 

Also,  the  corresponding  p.d.f.  is 


pf*)  ■ nW 


v42  + ff2/sinh2x 


(30) 


and  the  corresponding  variance  is 


2|  = 2L.  (£Othx 

9 'min  7 1 v 7 


sinh^x 


(31) 


Notice  that  many  values  of  x solve  equation  (28)  and  (29)  each 
one  corresponding  to  a local  minimum.  From  Figure  3,  it  is  clear 
that  °<plmin  decreases  as  x increases.  The  global  minimum  will 
correspond  to  the  highest  value  x that  solve  these  equations. 

5.  Remarks  on  the  Optimal  Solution 

(1)  Substituting  equation  (27)  in  equation  (2)  we  have 
for  the  system  equation 


d9  1+X2p 

3T  - A0  • H0AKfD^  - T+T^p 


a AKfj (9)  - HqKN*  (32) 
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where  N*  = Nf  + Nf  (1  -HQ)/H0 (1  + ^p) ; a*  = ill  ^ + t2 


fD(*)  ■ 

<P  +Ir  /sinh  x 


(33) 


and 


Aq  H0fVa 

fl(' : ) = Ate  ' Ctj+2‘  + A0/Aka)(a-1+H0T 


[-1  + M »42+"2/sinh2x] . 


(34) 


The  model  of  the  system  that  corresponds  to  these  equations  is  given 

in  Figure  4,  where  we  notice  that  channel  1 is  controlled  by  an 

odd  phase  detector  function  fp(4>),  while  channel  2 is  controlled 

by  an  even  phase  detector  function  f ^ (<)> ) . For  small  <j> , we  have 

fp(<f>)  - a^  where  a^  = sinh  x/n  « and  f^(40  - const.  That  is 

channel  2 does  not  respond  to  small  changes  in  <|>  and  the  system 

behaves  as  a first  order  PLL.  Noticing  that  the  phase  detector 

slope  a^  is  very  large  for  reasonable  value  of  x,  our  results 

resemble  similar  ones  obtained  by  Stiffler  [4]  and  Shaft  (5)  , 

namely,  the  phase  detector  function  that  causes  a first  order  PLL 

2 

system  to  have  minimum  is  sgn(sin  o) . However,  with  large 

or  PLL  unlocked,  the  system  acts  as  a second  order  PLL. 

(2)  Clearly  the  optimal  solution  depends  on  the  ratio 

u/B  (see  equation  (29)  which  itself  depends  on  Sf  (0)  (see 

*1 

equation  (9  )).  It  is  reasonable  to  assume  that  the  ratio  of 

the  noise  bandwidths  of  the  stochastic  processes  <J>  (t)  and 

f ^ (<J> ) - C4» ) is  proportional  to  the  ratio  of  their  variances 

2 2 

o^  and  o^  . In  this  case,  one  can  show  that  S^CO)  is 
1 2 

approximately  given  by  nNQ/2A  , where  n is  the  proportionality 
constant.  (Notice  that  for  the  regular  PLL  (i.e.  f^($)  ■ f(t)  ■ sin  ♦) 
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Viterbi  [2]  and  Lindsey  [1]  estimated  S^(0) 
this  estimate  for  S^.(0)  we  have 


by 


N°  ) 
tk:-> 


Using 


6/a 


N0 

l,0  4-tjA 


7 • 


(35) 


That  is  6/a  depends  on  the  noise  to  signal  ratio  of  the  input. 

(3)  Notice  that  the  optimal  solution  (see  equations  (28) 
and  (29))  depends  on  the  parameter  T^.  Recalling  that  this  para- 
meter (see  equation  (25))  must  vanish  for  the  standard  PLL.  To 
show  the  effect  of  T^  on  the  locking  time;  we  use  the  Fourier 
series  expansion  of  fj(4>)  in  equation  (32)  to  get 


d<j> 

at  = 


[a0  * 


(l+tip) a AKT, 

-Tprl  - h0ak£dw> 


Trrrp- 


(l+t*p)a 
-TTT—  AK 


HQKN* 


(36) 


where  we  also  used  the  fact  that  the  first  term  of  this  series 


is 


1/2 


' I b 


(<j))d<j>  * -T./27I,  and  gn  = £ a cos  n<J> . This 

1 u n=l  n 


equation  represents  a system  similar  to  that  of  Figure  4 but  with 
the  two  phase  detector  functions  f D (4> ) and  gD(<|>)  for  channel  1 
and  channel  2 respectively,  and  equivalent  input  signal  with 
phase  modulation 


1 U*T$P) 

Vp)  p lAo  + -TTxTF- 


AKaTj/2^) . 


(37) 


For  t >>  T^t  we  have  for  the  equivalent  frequency  input 
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d0  (t) 
eq  v J 

3t 


= V1 


aT^/2n 

Tjnr)- 


(38) 


Since  the  locking  time  of  a PLL  depends  on  the  incoming  signal 
frequency,  obviously  affects  this  time. 


6.  Computational  Results  and  Simulation 

As  it  was  shown  in  the  previous  section,  the  two  optimal 
phase  detector  functions  f(<J>)  and  f^(9)  depend  on  the  solution 
of  the  nonlinear  equations  (28)  and  (29).  Except  for  Tj  = 0, 
an  analytic  solution  is  difficult  to  obtain  and  hence  a numerical 
solution  is  obtained. 

In  fact  for  = 0,  it  is  possibleto  show  that  the  desired 
solution  to  these  nonlinear  equations  is 


a - 2(1-H0) 


(39) 


• and 


CA0AK)2a 

x “ nrnyv 


X >>  1. 


(40) 


In  this  case,  we  also  have 


o 


♦ ' min 


7,2  Cl-H0)H0B/a 

X (a0/ak)2 


(41) 


Notice  that  of  I . 

9 'min 

seems  unreasonable. 


decreases  as  AQ/AK  increases,  a result  which 
However,  using  the  loop  equation  (32),  we  can 


have 
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Aeff  4(1-H0)H03/u 

tit  = ' 


(42) 


That  is,  A rr/AK  decreases  as  A-/AK  increases,  and  the 
err  u 

2 

dependence  of  on  A^/AK  is  justified. 

With  the  values  of  a and  x from  equations  (40)  and  (41) 

the  optimal  phase  detector  functions  of  equations  (33)  and  (34)  are 

presented  in  Figures  5 and  6.  The  corresponding  optimal  p.d.f 

of  equation  (30)  is  given  in  Figure  7. 

For  T^  / 0 computer  search  can  be  used  in  solving 

equations  (28)  and  (29).  As  a numerical  example  consider  in 

Figure  8,  °^lm£n  as  a function  of  a/3,  (which  depends  on  S/N, 

equation  35)  with  system  parameters  r = 2,  Hg  = 0.02,  Tj  = 0.03 

and  for  An/AK  = 0.12  and  0.08.  Notice  that  for  a fixed  a/8, 

2 A0 

o , decreases  as  tv  increases,  similar  to  what  we  had  for  T.  = 0. 

<P  AK  .1 

2 

In  Figure  9,  °^lmin  as  a function  of  Tj  is  presented  with  system 

parameters  a/ 8 = 3 , r = 2 , Hq  = 0.01  and  for  Aq/AK  = 0.12  and 

2 

0.08.  Notice  that  for  a fixed  a/6,o  decreases  as  Tj  increases. 

The  theoretical  results  obtained  in  the  previous  sections  assumed 

that  S,  (0)  is  practically  independent  of  f , (4> ) (see  Remark  2 
tl 

of  Section  5).  If  this  assumption  is  true,  then  n should  be  in- 
dependent of  f^).  Thus,  to  justify  this  approximation, 

Figure  9 includes  simulation  results  of  as  a function  of  T^, 

with  Aq/AK  - 0.08  and  0.12,  and  for  n ■ 0.08.  Although  this 
value  of  n was  obtained  for  a single  point  on  the  curve,  the  small 
difference  along  the  curve  between  this  simulation  and  the  previously 
computed  results  do  support  the  assumption. 

In  Figure  10,  simulation  results  of  the  locking  time  as  a 
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function  of  T j in  a noise  free  system  are  presented.  The 

simulation  performed  for  optimal  system  whose  parameters  are 

A0 

r = 2,  H0  = 0.01  ^ = 0.12  and  u/p  = 3.  Locking  time  was 

i (1$  i ■ 4 

defined  to  be  the  time  for  the  system  to  reach  1^1  < 10 

Notice  that  the  locking  time  has  minimum  value.  Therefore  the 

design  goal  could  be  to  select  Tj  to  have  the  largest  value 

which  still  satisfies,  according  to  Figure  10,  the  locking  time 

requirement  (i.c.,  locking  time  < t ) . It  should  be  emphasized 

that  with  the  regular  PLL,  T^  is  necessarily  zero  and  there  is 

no  extra  freedom  in  selecting  the  locking  time  for  a given  AQ/AK. 

Previously,  the  assumption  that  Aq  is  known  was  used. 

This  is  not  true,  in  practice.  Hence,  it  is  important  to  discuss 

the  behaviour  of  a system,  optimally  designed  to  track  Ag  > when 

2 

it  is  actually  tracking  AQ  t Ag.  That  is,  to  find  of  a system 

having  phase  detector  functions  f|(<j>)  and  f*  (<)>)•  It  is 
practically  impossible  to  establish  analytically  the  p.d.f.  p(<f>) 

J J ft 

of  this  case.  Instead,  o^/o^  as  a function  of  Ag/Ag  was  calcuIated 

numerically.  The  results  are  depicted  in  Figures  11  and  12,  from 

2 2* 

which  it  is  noted  that  c^/o^  is  not  symmetric  with  respect  to 
AQ/Ag  =1.  It  increases  for  Ag/Ag  < 1 much  faster  than  for 
Aq/Aq  <1.  In  view  of  this,  we  may  conclude  that  for  Ag/AK  not 
known  a priori,  it  is  reasonable  to  design  the  system  for  the 
maximal  expected  incoming  frequency  error. 

7.  Practical  Implementation  of  the  Optimal  Phase  Detector 

Instead  of  the  optimal  forms  of  the  functions  given  on 
equations  (33)  and  (34)  or  Figures  5 and  6,  it  is  suggested  that 
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the  first  term  of  their  Fourier  expansions  be  used,  i.e., 

an 

f ^ (<l> ) - bjSin  <t>  and  f ^ ) ~ + a-^cos  4> , where  b^,  aQ  and  a ^ 

are  the  Fourier  coefficients.  To  examine  this  approximation, 

2 

Figure  13  compares  the  dependence  of  o^  on  Tj  of  this  system 
with  that  of  the  optimal  one.  It  is  noted  that  with  this  sub-optimal 
system,  the  variance  of  the  phase  error  is  not  much  larger  than 
that  of  the  optimal  one,  wfiich  justifies  the  usefulness  of  using 
these  easy  to  implement  phase  detector  functions.  Figure  14 
presents  a block  diagram  of  such  implementation. 
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Figure  5. 
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Model  of  the  PLL  with  two  phase  detectors 

Model  of  standard  PLL 
2 
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4 ' min 

Simplified  model  of  second  order  PLL  with  two  phase 
detector  functions 
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Figure  7. 

Figure  8. 
Figure  9. 

Figure  10. 
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The  optimal  probability  density  function  p(4)  with 


of  I as  a function  of  a/8  (SNR),  T.  f 0 

9 ' min  ’ 1 

2 

Simulation  and  theoretical  results  of  I min 
as  a function  of  Tj 

Locking  time  as  a function  of  Tj  in  a noise  free  system, 
using  simulation 


2* 
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as  a funct-on  of 
as  a function  of 
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* 
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Comparing  as  a function  of  T^  for  the  optimal 

system  with  a system  whose  phase  detector  functions 
are  fp(4)  = bjsin  4,  f ^ C4> ) = aQ  + ajcos  4 


Figure  14. 


The  sub-optimal  system 
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